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Abstract. A Golomb ruler is a sequence of distinct integers (tlie markings of the ruler) whose pairwise 
differences are distinct. Golomb rulers can be traced back to additive number theory in the 1930s and have 
attracted recent research activities on existence problems, such as the search for optimal Golomb rulers (those 
of minimal length given a fixed number of markings). Our goal is to enumerate Golomb rulers in a systematic 
way: we study 

gm{t) := # |x e : = xo < XI < • ■ • < x,„_i < x,„ = f, all Xj —xi^ distinctj , 

the number of Golomb rulers with m + 1 markings and length /. Our main result is that gm(t) is a quasipolyno- 
mial in t which satisfies a combinatorial reciprocity theorem: (— l)"'^'g,„(— ?) equals the number of rulers x of 
length t with m + 1 markings, each counted with its Golomb multiplicity, which measures how many combina- 
torially different Golomb rulers are in a small neighborhood of x. Our reciprocity theorem can be interpreted 
in terms of certain mixed graphs associated to Golomb rulers; in this language, it is reminiscent of Stanley's 
reciprocity theorem for chromatic polynomials. Thus in the second part of the paper we develop an analogue 
of Stanley's theorem to mixed graphs, which connects their chromatic polynomials to acyclic orientations. 



1. INTRODUCTION 

A Golomb ruler is a sequence of n distinct integers whose pairwise differences are distinct: one can 
picture an actual ruler with n markings having the property that all possible measurements are of distinct 
length. Golomb rulers have natural appUcations to phased array radio antennas (see, e.g., |2J, possibly 
the earliest reference on Golomb rulers, though not under that name), x-ray analysis of crystal structures 
(see, e.g., |5|), and error-correcting codes (see, e.g., |fT3l ). Mathematically they can be traced back to 
additive number theory in the 1930s ||9j [T3. In the more recent past, researchers have typically stud- 
ied existence problems, such as the search for optimal Golomb rulers (those of minimal length given a 
fixed number of markings), often with an eye toward computational complexity. Figure [T] shows an op- 
timal ruler of length 6. See, e.g., lU ISl [71 [Ml [161 and the parallel-search project on Golomb rulers at 
[iittp : / / www . dis tributed . net/Pro j ects. 

Our goal is to enumerate Golomb rulers in a systematic way. We define 

gm{t) := #{x € Z'"^^ : = ;co < xi < • • • < x^-i <Xm = t, all Xj —x^ distinct} , 

the number of Golomb rulers with m + I markings and length t. Most often it will be more convenient for 
us to express this counting function in the equivalent form 



(1) ^^(0=#<!zGZ^o: 



Zl +Z2H \-Zm =t, 

ZjeuZj + IjGvZ; for all dpcs U y d [m] 

where [m] := {1,2, . . . and dpcs stands for disjoint proper consecutive subsets, i.e., two sets of the form 
{a,a+ l,a + 2,. .. ,b}, {c,c+\,c + 2,. .. ,d} for some \ <a<b <c <d <m. 
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Figure 1. An optimal Golomb ruler with four markings. 



Motivated by ([T]), we call a point z G Z>q with zi+Z2-\ hZm = ? a ruler of length t; thus a Golomb 

ruler is a ruler with positive entries satisfying Y^jeu^j / L;ey for all dpcs ?7, V C [m]. 

We define a rea/ Golomb ruler (with m + 1 markings and length t) as a vector z G M>q satisfying 

Zi + Z2 H VZm = t and XI 7^ £ Zy for all dpcs U,V (Z[m\. 

JeU jev 

Two real Golomb rulers z, w € M'"q are combinatorially equivalent if for any dpcs U ,V C [m], 

jeu jev jeu jeV 

in plain English: if their possible measurements satisfy the same order relations. We define the Golomb 
multiplicity of a ruler z e Z>q to be the number of combinatorially different real Golomb rulers in an e- 
neighborhood of z (viewed as a point in W"), for sufficiently small £ > 0. Thus a Golomb ruler is a ruler 
with Golomb multiplicity 1 . 

A quasipolynomial of degree <i is a function of the form Cd{t)t'^ + Cd-i{t)t'^^^ H \-co{t), where the 

Cj{t) are periodic functions in t. The 1cm of their periods is called the period of this quasipolynomial. Our 
main result about gm{t) is the following combinatorial reciprocity theorem. 

Theorem 1. The Golomb counting function gm{t) is a quasipolynomial in t of degree m—\ with leading 
coefficient -j^^^riy- Its evaluation ( — l)'"^'gm(— equals the number of rulers in Z>q of length t, each 

counted with its Golomb multiplicity. Furthermore, ( — l)'"^^gm(0) equals the number of combinatorially 
different Golomb rulers with m+\ markings. 

Example 1. Let m = 3. Using James B. Shearer's Fortran code for Golomb ruler^ we computed the values 
of g3 (f ) for 6 < ? < 35. Note that if z is a Golomb ruler with m + 1 markings, then there is a complementary 
Golomb ruler z' given by z- = Zm-i, and (z')' = z. The rulers z and z' are always distinct unless m = \. 
Shearer's code computes only one of each complementary pair, but in Table [T] we double the output values 
to account for both. We also note that by inspection, g3{t) = for ? < 6. 

In Section |2] we show that the period of the quasipolynomial function gi{t) divides 12 and that its leading 
term is ^t^. Thus from the 24 values g3(0),g3(l), . . . ,g^{12>) we can obtain all of the coefficients by interpo- 
lation. The remaining 12 values g3(24), . . . ,g3(35) are consistent with the result of this computation, which 
is: 



Available at http: //www. research, ibm. com/people/s/ shearer/programs/grsl . txt 
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Table 1. The Golomb counting function gT,{t) for 6 < f < 35. 
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In particular, the coefficient c\{t) has period 2. The coefficient C()(f) has period 12 but obeys the same 
formula for t = j (mod 12) as for f = —j (mod 12) for each / Note also that ^3(0) = 10; Theorem [T] 
predicts that there are ten combinatorially different Golomb rulers, as we will see in Section |2] 

We prove Theorem [T] geometrically, making use of the machinery of inside-out poly topes 0- This ap- 
proach leads us to associate a mixed graph (i.e., a graph that may contain both undirected and directed edges; 
see, e.g., [12]) to the set of all Golomb rulers with a fixed number of markings. In this language. Theorem 
[T] has an interpretation in terms of acyclic orientations of mixed graphs, which might be of independent 
interest (Theorem |4] and Corollary [5l proved in Section |3]). This reinterpretation is reminiscent of Stanley's 
reciprocity theorem for the chromatic polynomial of a graph fT9l. This leads us to a natural analogue of 
Stanley's theorem for (general) mixed graphs, as follows. 

We write a mixed graph G as G = {V,E,A) where E represents the undirected edges and A the directed 
edges0 A t-coloring of a mixed graph G is a map c : V — )• [f]. Such a f-coloring is proper if 

• c(v) / c{w) for all {v,^} G E, and 

• c(v) < c{w) for all (v, w) € A. 

As with undirected graphs, the chromatic number is the minimum f such that G can be f-colored. This 
parameter was introduced for mixed graphs by Hansen, Kuplinsky, and de Werra ifTTl . who discuss bounds, 
algorithms, and applications to scheduling problems. Let Xcit) denote the number of proper f -colorings of 
G. Sotskov, Tanaev, and Werner [.18.1 showed that this function (if not identically zero) is a polynomial of 
degree \V\ and computed the two leading coefficients. 

An orientation of a mixed graph is obtained by keeping each directed edge and orienting each undirected 
edge; an orientation is acyclic if it does not contain any coherently-oriented cycles. A coloring c and an 
orientation of G are compatible if 

c{u) < c(v) <^=^ M — >■ V. 
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We assume that G is simple: that is, that it contains neither muhiple edges nor loops. 
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Note that c may not be a proper coloring. 

Here is our generalization of Stanley's chromatic reciprocity theorem to mixed graphs. Our approach (in 
SectionlU also yields a new proof of polynomiality of Zcit)- 

Theorem 2. Let G be a mixed graph. For a non-negative integer t, (— l)'^'j^^G(~0 ^(luals the number of 
t-colorings of G, each counted with multiplicity equal to the number of compatible acyclic orientations of 
G. 

Corollary 3. For G a mixed graph, (— 1)'^'^g(~1) equals the number of acyclic orientations of G. 

In Section[5]we end with some open questions about Golomb rulers suggested by our geometric approach. 

2. INSIDE-OUT POLYTOPES 

It is a short step to interpret a Golomb ruler in its measurement representation (as in ([T])) as an integer 
lattice point confined to the positive orthant in and the affine space defined through 

Z\ +Z2^ VZ,n=t. 

More precisely, a Golomb ruler with m+\ markings and length f is a lattice point in the f'th dilate of 

A°,:={zGR™o: Zi+Z2 + --- + Zm = l}, 

the (m — 1) -dimensional open standard simplex Uving in M'". However, this description does not yet capture 
the distinctness condition: a Golomb-ruler lattice point must also avoid the hyperplanes in W" given by 
equations of the form 

yec/ jev 

for all dpcs U C [m]. Let '^^ be the collection of all such hyperplanes. Thus we have built a geometric 
setting in which to compute the Golomb counting function as 

(2) g,„(0 = #(?A° \$f„,nz'"). 

In the language of [4], (A„,,?^„,) is an inside-out polytop^ and gm{t) is the (open) Ehrhart quasipolynomial 
of (Am,§f„j). Figure |2] shows this geometric setting (viewed in the plane Zi+Z2 + Z3 = 1) for the case m = 3. 




Zl =Z2 Zi =Z2+Z3 



Figure 2. The inside-out poly tope (As,'^!,). 



^To be 100% precise, f4\ would compute g,„ [t) as # (A°, \ Sf,„ n |Z'") , but this is equivalent as the hyperplanes in are linear. 
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The vertices of this inside-out polytope are 

{OAl),i^A^) ,{0,^,^) ,{1,0,0), {^,^,0) , (0,1,0). 

The 1cm of the denominators among all of these coordinates is twelve. This shows that the period of the 
quasipolynomial g3{t) divides twelve, as claimed in Example [U We can also see that dissects A3 into ten 
polygons, which correspond to the combinatorially different Golomb rulers. 

The general setup of an inside-out polytope {£P,J^) consists of a rational polytope ^ and a rational 
hyperplane arrangement J^f in M^; that is, the linear equations and inequalities defining ^ and have 
integer coefficients. The goal is to compute the counting function 

L°^Mt) :=#((^\jr)niz''), 

and it follows from Ehrhart's theory of counting lattice points in dilates of rational polytopes |l3l. (81 that 
this function is a quasipolynomial in t whose degree is dim{P), whose (constant) leading coefficient is the 
normalized lattice volume of P, and whose period divides the 1cm of all denominators that appear in the 
coordinates of the vertices of {^,J^). Furthermore, [4J established the reciprocity theorem 

(3) L°g^oj^{—t) = {-\)'^^'^^L-^j^{t), 

where and J!P denote the interior and closure of respectively, and 

(4) L^^,^{t) := 52 muk^i2.^^(m) 

where mult;^ ^^(m) denotes the number of closed regions of {l^,J^) containing m. (A region of J^) 
is a connected component of \ J^f; a closed region is the closure of a region.) It follows from Ehrhart's 
work and (HJl that 

(5) L^,,^(0) =# regions of 

See PI for this and several more properties of inside-out polytopes. The concept of inside-out polytopes has 
been applied to a number of combinatorial settings; at the heart of any such application is an interpretation 
of the regions of J^), which we will now give for the Golomb inside-out polytope (Am,^^,,,)- 

Proof of TheoremUl The first statement follows immediately from the fact that {Am,'^m) is a rational inside- 
out polytope and A„, is a unimodular {m — 1) -dimensional simplex. 

Viewing ^ from Q, we know that (— l)'"^^g,„(— f) equals the number of rulers in An jZ'", each counted 
with multiplicity given by the number of closed regions of (A„,,?^„j) it lies in. These regions, in turn, are 
defined by inequalities of the form 

jeu jev 

for some dpcs U,V C [m], and thus the multiplicity of a ruler z is precisely given by the number of combi- 
natorially different real Golomb rulers in a neighborhood of z. 

The last statement of Theorem[T]follows from ([5]). □ 

3. The Regions of (A^,^™) 

For each positive integer m, we define the Golomb graph F,,, to be a mixed graph whose vertices are all 
proper consecutive subsets of [m]. The underlying graph is complete and an edge UV is directed {U V) if 
and only ifU CV. All other edges are undirected. Acyclic orientations of a Golomb graph allow us to give 
an interpretation of the regions of (A„,,^„) in the following sense. 
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Theorem 4. The regions of the Golomb inside-out poly tope {Am,^m) fl^re in one-to-one correspondence with 
the acyclic orientations of the Golomb graph F,,, that satisfy the relation 

(6) A^B U->V. 

for all proper consecutive subsets A and B of [m] of the form A = UUW and B = VU W for some nonempty 
disjoint sets U,V,W 




Figure 3. The mixed graph F3 (left) and an acycUc orientation of F3 (right). 



Proof. Let Rhea region of ( A,„ , '^,„ ) . Thus R is the intersection of A„, with half spaces defined by inequalities 
of the form 

(7) L^j<L 

jeu jev 

for some dpcs U,V C [m]. So UV is an undirected edge of F,„; assign it the orientation U ^V. 

The directed edges of F^ are of the form UV for U C V, so we still need to orient edges of the form 
AB where A = U UW and B = V UW for some nonempty disjoint sets U,V,W. We orient these according 
to (O (note that the corresponding half space defined through Y^jeA^j < HjeB^j is identical to that given by 
([7])). The resulting orientation of F^ is acyclic since a coherently-oriented cycle gives rise to the nonsensical 
inequality Y^jeM^j < HjemZj for some multiset M. 

Conversely, suppose we are given an acyclic orientation of F^ that satisfies Index the variables yu 
of a vector in M"'('"+i)/2^i by the proper consecutive subsets U of [m] and consider the braid arrangement 
given by the hyperplanes yu =yv in 0/2-1 ^ famous observation of Greene [10 ] that the acyclic 

orientations of a graph are in one-to-one correspondence with the region of the corresponding graphical 
arrangement (in our case, the above braid arrangement, since our underlying graph is complete). Now 
consider the function M'»('«+i)/2-i ^ m.'" defined through 

This is a linear map of full rank, i.e., the regions of the braid arrangement, being full-dimensional sets in the 
domain, are mapped to full-dimensional sets in the range. Furthermore, since the all-ones vector is in each 
closed region, we can project further to the affine space defined by zi + Z2 H \-Zm = f- D 

This theorem allows us to rephrase Theorem [Din terms of acyclic orientations of F^. We call a ruler 
z G Z>Q and an orientation of F„, compatible if 



jeu jev 



U ^V. 
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Corollary 5. The evaluation ( — equals the number of rulers in Z'"q of length t, each counted 
with multiplicity equal to the number of compatible acyclic orientations of that satisfy Further- 
more, ( — l)'"^'gm(0), the number of combinatorially different Golomb rulers, equals the number of acyclic 
orientations ofr,„ that satisfy (O. 

4. Mixed Graph Coloring Reciprocity 

Our proof of Theorem[2]closely follows the proof of Stanley's chromatic reciprocity theorem given in f?*]. 
Given a mixed graph G = {V,E,A) with V = [n], define the polytope 

^(G) := {x G [0, 1]" : xj < Xk whenever {j,k) G A] 

and the hyperplane arrangement J^{G) consisting of the hyperplanes xj = Xk for all {vj,Vk} G E. The 
following result is the mixed-graph analogue of ifTOl . 

Proposition 6. The regions ofM"\J^{G) that intersect ^{G) are in bijection with the acyclic orientations 
ofG. 

Proof. Let R be an open region of J^(G) that intersects ^(G), hence also intersects ,0^{G)°. Then for each 
hyperplane x,- = Xj in .j^{G), R is entirely on one side, hence determines an orientation a as follows: for 
each edge {/, 7} G E, choose / — > j if x,- < xj in R, or j — )• / if x,- > Xj in R. In the set R n l3^{G)° , we also 
have Xi < xj for all (/, j) G A. Since R, being open, intersects the interior of ^{G)°, all of these inequalities 
must be consistent. That is, a is acylic. 

Since two different regions must be on opposite sides of at least one hypeiplane, they determine different 
orientations. Finally, if a is an ayclic orientation of G, choose any total order n of the vertices that is 
consistent with a, including the fixed directions on the arrows A. Then if R is the region containing the 
points of M" that satisfy x^^j < • • • < X;^,, > then R maps to a under the rule given above. That is, the map from 
regions to orientations is both injective and surjective. □ 

We say that G is acylic if the directed graph (G,A) contains no directed cycles. Note that in this case there 
exists at least one total order on the vertices compatible with the directed edges. Orienting the undirected 
edges according to such an order, we can always obtain an acyclic orientation on all of G. Note that if G 
is not acyclic, then there exist no proper ^-colorings for any t. On the other hand, if G is acyclic then there 
exist proper ^-colorings for all sufficiently large t, and in particular for all t >n. 

Lemma 7. If the mixed graph G is acyclic then the polytope ^{G) is full dimensional. 

Proof. The polytope ^{G) depends only on the directed graph D = {V,A). Since D is acyclic, there is a 
total order ii < ii < ■ ■ • < in on [n\ such that j <k for all {vj^Vk) G A. Thus 

^{G) 2 {x G M" : < X;, < x,^ < • • • < x,„ < 1} 

which is a full-dimensional simplex in M", so .!^{G) is full-dimensional. □ 

Lemma 8. For any acyclic mixed graph G, the inside-out polytope (G) , Jif (G)) is integral. 

Proof. Any closed region of {j^{G),J^{G)) is defined by a collection of inequahties of the formxy — x<; < 
and Xj > 0. In particular, the left-hand sides of these inequalities are all vectors of the form ej — tk or — ey. 
Any matrix whose rows are all such vectors is totally unimodular. This implies that all vertices of the region 
are integral ifTSl . □ 

Proof of Theorem^ Set ^ := ^{G) and M' := M'{G) . Using the identification of maps c -.V 
with lattice points in M", we see that for f > 1, c is a — 1) -coloring (respecting the directed edges A) if 
c G t^°{G) and c ^J^. That is, Xcit - 1) = which equals (-l)"£^_jr(-0 by ©. Using 

LemmalU we see that the functions E^ ,^^ and E°^o ^ are polynomials and hence the reciprocity extends to 
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all values of t. In particular, for s a positive integer we have Xg{—s — 1) = E°^o j^{—s) = {—\)"'E^^j^{s). 
Now Eg>^^{s) counts each point c G si^ with multiplicity equal to the number of closed regions of that 
intersect the interior of ^ and contain c. By Proposition [6l regions correspond to acyclic orientations of 
G. The closed region corresponding to an orientation a contains c if and only if c{u) < c(v) for each edge 
M ^ V in a. □ 

Example 2. Let G = ([3], {{1,3}, {2, 3}}, {(1,2)}), a triangle with one directed edge and two undirected 
edges. Of the four orientations of G, three are acyclic and each of these three determines a total order -< on 
the vertices. These three orders are: 

1^2^3,1^3^2,3^1^2. 
The polytope ^(G) is the half-cube given by the five inequalities 

0<Xi <X2 < 1,0<X3 < 1. 

The arrangement Jif{G) consists of the two hyperplanes xi = X3 and X2 = X3, and the inside-out polytope 
(^(G), J^(G)) has three open regions in bijection with the acyclic orientations of G: 

^123 = {x G : < xi < X2 < X3 < 1 } 
^132 = {x e : <xi <X3 <X2 < 1} 
^312 = {x G : <X3 <xi <X2 < 1} . 

For this graph it is easy to compute Xg{^) directly. In any coloring c, we must choose c(l) < c(2) and 
then we can choose c(3) to be any remaining color Thus there are (2)(f — 2) = '^'^^}^'^^^ f-colorings for 

any non-negative integer t. Then (— l)^j^G(~l) = — = 3, which is the number of regions of 

(=^(G),^(G)) as stated in TheoremE] 

As we know from the proof of Theorem [2l we can also use Ehrhart reciprocity to interpret the evaluation 
of Xg at other negative integers. Specifically, —Xci—t) is the number of lattice points in {t — l)^ with 
multiplicities, where the multiplicity of each point is the number of closed regions of the inside-out polytope 
that contain it. For instance, -Xg{-2) = (-1)^^^^^^^^ = 12 and -Xg{-^) = (-1)^^^^^^^^ = 30. 
In Table |2j we list the lattice points in ^ and 2^ and the regions that contain each point. Indeed, the 
multiplicities add up to 12 and 30, respectively. 

Note that in this example, J^{G) is not transverse to J^(G). The two hyperplanes xi = X3,X2 = X3 in 
J^{G) meet in a flat given by xi = X2 = X3. This flat intersects only the boundary of ^ since the equation 
xi = X2 is a supporting hyperplane of 



5. Open Problems 

The search for optimal Golomb rulers remains a challenging open problem. This problem is related to 
our paper: the smallest positive integer that is not a root of g,„ (t) gives the length of an optimal Golomb ruler 
with m + I markings. Our methods present new open problems, from the computation of g,„{t) to specific 
questions, e.g., about the period of g,„{t) and the number of combinatorially different Golomb rulers with 
a given number of markings (this sequence starts with 1, 2, 10, 114, 2608, and 107498). Our reciprocity 
theorem for chromatic polynomials of mixed graphs suggest open problems about these polynomials, from 
classification problems (which polynomials are chromatic polynomials of mixed graphs?) to interpretations 
of and relations between its coefficients. 
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point z G ^ n 


regions containing z 


multiplicity of z 


(0,0,0) 


^I23,-'?I32,-^3I2 


3 


(0,0,1) 


■^^123 


1 


(0,1,0) 


^132,^312 


2 


(0,1,1) 


^123,^132 


2 


(1,1,0) 


^312 


1 


(1,1,1) 


■'?I23,-'?I32,^3I2 


3 




point z G 2^nZ^ 


regions containing z 


multiplicity of z 


(0,0,0) 


^123,^132,^312 


3 


(0,0,1) 


^123 


1 


(0,0,2) 


^123 


1 


(0,1,0) 


^132,^312 


2 


(0,1,1) 




2 


(0,1,2) 


Rm 


1 


(0,2,0) 


^132,^312 


2 


(0,2,1) 


^132 


1 


(0,2,2) 


■/?132,-'?I32 


2 


(1,1,0) 


R312 


1 


(1,1,1) 


^123,^132,^312 


3 


(1,1,2) 


/?123 


1 


(1,2,0) 


■'?3I2 


1 


(1,2,1) 


^132,^312 


2 


(1,2,2) 


■'?I23,^I32 


2 


(2,2,0) 


R312 


1 


(2,2,1) 


R312 


1 


(2,2,2) 


^123,^132,^312 


3 



Table 2. Lattice points in and 2^ and the regions that contain each point. 
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